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Abstract 

An exactly integrable symplectic correspondence is derived which in a continuum limit 



Qh. leads to the equations of motion of the relativistic generalization of the Calogero-Moser 

' system, that was introduced for the first time by Ruijsenaars and Schneider. For the 

^ . discrete-time model the equations of motion take the form of Bethe Ansatz equations for 

^ I the inhomogeneous spin-^ Heisenberg magnet. We present a Lax pair, the symplectic 



structure and prove the involutivity of the invariants. Exact solutions are investigated in 
the rational and hyperbolic (trigonometric) limits of the system that is given in terms of 
elliptic functions. These solutions are connected with discrete soliton equations. The re- 
sults obtained allow us to consider the Bethe Ansatz equations as ones giving an integrable 
symplectic correspondence mixing the parameters of the quantum integrable system and 
the parameters of the corresponding Bethe wavefunction. 
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1 Introduction 



In some previous papers, |||, ^, cf. also an exact time-discretization of the famous 
Calogero-Moser (CM) model, was introduced and investigated. The discrete model 

is an integrable symplectic correspondence, (for a definition cf. that in a well-defined 
continuum limit yields the classical equations of motion of the CM system^. A few years 
ago Ruijsenaars and Schneider introduced in cf. also ||ll|, a relativistic variant of 
the CM model, which is a parameter-deformation of the original model. The equations of 
motion of this system in its generic (elliptic) form read 

qi = qiqjv{qi - qj) , i = l,...,N , (1.1a) 
where the potential v{x) is given by 

"<^-' - wrh> • <"^' 

in which p(x) = p{x\uji,uj2) is the Weierstrass P-function, 2ci;i,2 being a pair of periods. 
(These are the equations of motion in the form given by Bruschi and Calogero in []l^]). 
This multi-particle model is also integrable, and carries a representation of the Poincare 
algebra in two dimensions. Moreover, a large number of the characteristics of the CM 
model are generalized in a natural way to the relativistic case, such as the existence of a 
Lax pair, a sufficient number of integrals of the motion in involution, and exact solution 
schemes in special cases. 

In view of the results in [|l|-[^], a natural question to ask is whether there exists a 
time-discrete version of the model ([1.1|). We have to remark that on the quantum level 
the transition from the usual CM system to its relativistic counterpart, already amounts 
to a discretization (or g-deformation). In fact, the relevant operators in the quantum 
relativistic model are commuting difference operators, rather than differential operators, 



cf. |14, 15|. In view of the importance of quantum Calogero-Moser type of models in 

the context of representation theory, and in particular in connection with the Knizhnik- 

Zamolodchikov (KZ) equations as has been revealed by a large amount of recent work, cf. 

^It should be noted that the discrete CM model can be inferred also from the Backlund transformations 
for the continuous CM model, that were presented in M. 



e.g. [16|-|l^, these difference operators yield new interesting connections with g-special 



functions. 



In this paper we will introduce a discrete-time version of the model (1.1), in the 



form of an integrable symplectic correspondence (i.e. multi-valued map) which goes in a 
continuum limit to the original model. This amounts to one more parameter-deformation 
of the CM model: apart from the "spatial" discretisation (encoded in the parameter 
A), the discretisation of time constitutes another deformation where the finite step-size 
in time enters as the new parameter. The construction is based on an Ansatz for a 
Lax pair, together with an elliptic version of the Lagrange interpolation formula. We 
will demonstrate the integr ability of this mapping along the lines of ref. |l^, i.e. 
integrability in the sense of Liouville. In contrast to the previous paper, Q, where the 
particle model is related to pole solutions of discrete soliton equations, in particular of 
the lattice Kadomtsev-Petviashvili (KP) equation, the relativistic case is related to soliton 
solutions of this type of discrete equations. Finally, we point out the intriguing resemblance 
between the equations of motion of the discrete particle model and Bethe Ansatz equations 
for integrable spin-^ quantum chains of the XYZ Heisenberg model (see ||3^, |39||). 

2 Derivation of the Model 

We will start from an Ansatz for a Lax pair of the form 

N 

L^= hihj^,^{qi - Qj + X)eij , (2.1a) 
Mk = Y1 hihj^^iqi - qj + X)eij . (2.1b) 



In (2.1) the qi denote the particle positions, and the hi are auxiliary variables which 
we specify later. The tilde is a shorthand notation for the discrete-time shift, i.e. for 
qi{n) = qi we write qi{n + 1) = q^, and qi{n — 1) = q^. The variable k is the additional 
spectral parameter, whereas A is a parameter of the system. The matrices Cij are the 
standard elementary matrices whose entries are given by {eij)k£ = Si,k^j/- The function 



$K is defined as 



. , (2.2) 



where cr(x) is the Weierstrass sigma-function, given by 

1 T „1 

(2.3) 



<t(x) = X TT (1 ) exp 



X 1 X , 9 

+ o( ) 

'■^ki ^ ^ki 



with uJki = 2kuJi + 2iuJ2 and 2u;i 2 being a fixed pair of the primitive periods, see e.g. [20|. 
The relations between the Weierstrass elhptic functions are given by 

C(x) = ^ , p(x) = -C'(x) , (2.4) 

where ct(x) and ^{x) are odd functions and p{x) is an even function of its argument. We 
recall also that the cr{x) is an entire function, and C{x) is a meromorphic function having 
simple poles at Uki, both being quasi-periodic, obeying 

Cix + 2001,2) = Cix) + 2??i,2 , a{x + 2a;i,2) = -a(x)e2''i.2(^+'^i,2) ^ 

in which r/1^2 satisfy r]ii02 — t]2u:i = ^ > whereas p{x) is periodic. From an algebraic 
point of view, the most important property of these elliptic functions is the existence of a 
number of functional relations, the most fundamental being 

C a +C /3 +C 7 -C « + /3 + 7 = ) ^ .J. ) ^ /' L X • 2.5 

From this relation, one can basically derive all important identities for the Weierstrass 
elliptic functions. For our purpose it can be recast into the form 

^>,(x)^>,(y) = $,(x + y) [C(k) + C(x) + C(y) -C{k + x + y)] . (2.6) 

The Ansatz for the Lax pair (^]^) is a natural one, in view of the fact that the matrix 
corresponds to the Lax matrix of the continuum Ruijsenaars- Schneider (RS) model, 

cf. [13 1, whereas the matrix is a natural choice by comparison with the earlier results 

obtained in Q, cf. also Q, for the discrete elliptic CM model. 

Let us now consider the compatibility of the system ( p.l| ) . Assuming the Lax equation 

L^M^ = M^L^ , (2.7) 

which implies the isospectrality of the discrete flow of the Lax matrix L^j we have from 
(^ that 

J2 + Ciqi-qi + A) + Ciqe- gj + A) - ((« + 2A + - qj)] 

= II ^£ [C('t) + C(gi - + A) + C(g^ - qj + A) - ((/^ + 2A + - qj)\ . 



Noting that the conservation law trL = trL imphes: 

E^' = E^'' (2-8) 
e. I. 

we have the identity 

E [^^C(% - + A) - hlC,{qi -qe + \) = - ^ [^IC(*^ - Qj + A) - /i|C(® - + A) , 

(2.9) 

for all i,j = 1,...,N: consequently, both sides of ( |2.9D must be independent of the 
(external) particle label. Thus, we find a coupled system of equations in terms of the 
variables hi,qi of the form 

E [^kiQi -qi-^)- hjCiQi -Qi- A)] = -p , (2.10a) 



e 

hj CiQi - qe +A) - hjCiQi -qg + X) 



E 



(2.10b) 



where p does not carry a particle label. In principle p can still depend on the discrete 
time-variable, but we will mostly consider it to be constant: a different choice will be 
taken only in Section 6. In that case, by eliminating the variables hi from ( |2.9| ), we get a 
closed set of equations in terms of the qi. 

In order to derive this closed set of equations for the variables qi we will make use of 
an elliptic version of the Lagrange interpolation formula, that was derived in [^. In fact, 
we have the following statement: 

Lemma: Consider 2N noncoinciding complex numbers Xi,yi (/ = 1, . . . , N). Then, the 
following formula holds true: 

/=i - yii f^i ni=i (^{yi - yk) 

k^l 

N 

where ^{yi - x;) = , 
(=1 

and where x stands for any one of the zeroes xi of the function on the left-hand side. 

Eq. ( |2.11| ) can be derived from an elliptic version of the Cauchy identity for the 
determinants of the matrix {^^{xi — yj))- This Cauchy identity was proven and used in 
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2[ to establish the commutativity of the quantum integrals for the continuum RS model. 



Actually, this identity goes back to Frobenius, |21], and can be proven also by purely 
combinatoric means starting from the fundamental identity ( |2.5D . A statement similar to 



the one in the lemma was formulated in |22|. 

What is essential in ( ^.11 ) is the independence of the choice of x, which can be easily 
demonstrated by using the identity 







(2.12) 



which can be found e.g. in |23], p. 451. 

Using now eq. ( p. Ill ) and applying the lemma to the elliptic fractional "polynomial" 



N 

n 

k=l 



cr{C -qk + A)cr(^ -Qk-X) 



<^{i - Qk)(j{i - Qk) 

and taking C = Qj ~ respectively x = qi + X , we can by comparing with eq. (^) make 
the identifications 

n^=i o-{q£ -qk + X)(y{qi -qk- X) 



hi 



-P 



Uk^e(^{qe- qk) Uk=i(^{q> 



qk) 



72 _ Y\k=i f^iqe -qk + >^)(y{qt -qk- X) 

— p 



Uk=ii(^{qi-qk) Uk=i(^{q 



(2.13a) 
(2.13b) 



qk) 



1,...,N 



from which we obtain the following system of equations 

N r ^^ N ^ili " ^fc) 0-{qi - Qk +A) 

P -TT crjqe - gfc + A) ^ -pj- ~ 

P criqi -qk- X) j^}^ a{qi - qk) cr{qi - qk - X) 



l,...,N . 



(2.14) 



Eqs. (|1|) can be considered to be a product version of (p.9|), and is a system of N 
equations for + 1 unknowns, qi,. . . ,qi\i and p. There is no equation for p separately, 
and thus it should be a priori given in order to get a closed set of equations. The most 
natural choice is the one for which p is constant, i.e. independent of the discrete time- 
variable. For convenience we will take it equal to the fixed value p = — cr(A)~^, but if 
we are interested in a continuum limit from ( 2.10| ) we should take it of the order of the 
reciprocal of the discrete-time step. As already mentioned, a dynamical choice of p arises 



at the end of Section 6 from the apphcation of ( 2.14| ) to the Bethe Ansatz equations for 
the XYZ Heisenberg chain. 



Thus, taking p/ p to be equal to unity in (2.14), we obtain the equations of motion 
of what we would like to call the discrete Ruijsenaars-Schneider model. It is given by 
a coupled set of algebraic equations, which, in fact, resemble closely the Bethe Ansatz 
equations (BAE's) for certain integrable quantum chains with impurities. We will make 
this connection more precise in Section 6. Here, we look at eqs. ( |2.14 ) as defining a 



discrete dynamical system, which amounts to an integrable symplectic correspondence, 
i.e. a multi- valued map in the sense of ref. Q. 

Before proving the symplecticity and integrability of the correspondence, let us show 



first that the interpretation of ( 2.14 ) as a discrete version of the RS model is justified, by 
showing that an appropriate continuum limit will give us back the equations (jLlD of the 
continuum model. In fact, taking 

. 1 2- 

gfe 1-^ gfc - A + egfc + -e gfc + . . . , 

1 2- 

£fc 1-^ gfc + A - egfc + -e gfc + . . . , 
and developing with respect to the small parameter e, we obtain for the leading term 

- = Y.ik [2C(g£ - qk) - C{qe -Qk + X)- C{qe - qk - A)] , (2.15) 

k^e 

from which we recover eqs. ( |l.lD using the relation 



C(..A,-C(.)-C(A)4i^:|^. (2.10) 

Let us finally remark that the non-relativistic limit is obtained by letting the parameter 
A go to zero. In fact, in this limit the variables hf will behave as 

h^^l + xlj2C{qi-qk)-Y.((^^-^k)] +0{X^) . (2.17) 

k=l J 

As A —> 0, we will obtain the matrix of the non-relativistic case immediately, whereas 
for we thus find 

^ X'^ + Y,Pie^i + Y,^Mi - qj^ij + C(A) , (2.18) 



where the pi are the expressions between brackets on the right-hand side of the arrow in 
( |2.17| ). Thus, we recover the Lax representation of the non-relativistic case as given in Q. 



3 Symplectic Structure and Integrability 

In order to demonstrate the integrabihty of the discrete-time RS model, given by the 



equations of motion (2.14) together with p = p , we need to establish an invariant 
symplectic structure. Following the philosophy of previous papers, (cf. e.g. and 
references therein), this can be assessed on the basis of an action principle, cf. also ||8|, |19|. 



It is easy to note that an action for the discrete equations ( 2.14| ) is given by 



S = Y.C{q{n),q{n+l)) , (3.1a) 
n 

in which 

TV TV 

^Q, q) = ~ ^fc) - file - 9fc - A)] - f{qi -Qk + X) , (3.1b) 



k,£=l k,e=i 



with the function f{x) given by 

fix) = I log\a{CM ■ (3.2) 



Let us point out that the function / appearing here can be considered to be an elliptic 
version of the Euler dilogarithm function. The discrete Euler-Lagrange equations 

1^ + 1^ = , i=l,...,N, (3.3) 
dqc dqe 



are easily seen to lead to eq. ( 2.14 ). The canonical momenta pc are found from 

dC ^ ^ 
P£ = ^ = y] (- log \cr{qe -qk)\+ log \cr{qe - Qfc + A)|) . (3.4) 

As a consequence we have that the symplectic form (7 = J2k ^Pk ^ dQk is invariant 
under the correspondence, which implies that any branch of the correspondence defines a 
canonical transformation with respect to the standard Poisson brackets given by 

{pk,q£} = Sk,i , {pk,Pi} = {qk,q£} = ■ (3.5) 
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We will now show that in terms of the canonical variables, the Lax matrix in ( ^.Ip 
has the same form as the one of the continuum RS model. In fact, expressing the variables 
hf in terms of the canonical variables, we obtain 

hj = e- n ^^^^ , (3.6) 



which we will refer to as the Ruijsenaars variables, |1C]. In terms of these variables we 
have the Poisson brackets 

{<ik, 9^} = , {log hi , qi} = 5k/ , 

{log hi , log hj} = C{qk-qe + A) + C{qk - qe - X) - 2C(% - qe) , k £ . (3.7) 

Now, there remains not so much work left to establish the integrability. In fact, having 
expressed the Lax matrix in terms of the original variables, and having verified that 
it takes the same form as in the continuous case, we can rely on the proof of involutivity 
that was produced by Ruijsenaars (in Appendix A of |[lO[| ) to assess the involutivity of the 
invariants of the discrete model as well. Naturally, the invariants of the correspondence, 
given by 

Ik = tTL^ = Ik , k = l,...,N, (3.8) 
as functions of the canonical variables are the same as in the continuum model. 
Thus, we obtain the statement of involutivity of the invariants, i.e., 

{h,k} = {Tr{L'),Ti{L')} = for all fe,£ = l,2, (3.9) 

It is a theorem by Veselov that the involutivity of the invariants lead to the linearization 
of the discrete flow on tori, similar to the continuous-time situation, cf. (For a simple 
argument see also |jl9|). 



4 Exact Solutions 



Let us now consider exact solutions of the discrete RS model in two special cases: i) the 
rational limit, ii) the hyperbolic (trigonometric) limit. Here, again, we restrict ourselves 
to the case p = p , which will then lead to an explicit integration scheme from the Lax 
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pair, as we shall demonstrate below. If, more generally, we do not fix p, the initial value 
problem becomes implicit, having to take into account initial values for the p also. We 
shall not pursue this latter case in this paper. 

The main line of reasoning in both limiting cases is very similar to the one followed in 
[0, ^ for the non-relativistic discrete CM model. In the discrete situation the initial value 
problem is posed by asking to find {qi{n)} for given initial data {(?j(0)} and {(?j(l) = qi{0)}, 
where qi{n) denotes the position of the particles Qi at the nth time-step. The solution to 
this problem is obtained for both the rational as well as the hyperbolic (trigonometric) 
limits by solving a secular problem, namely by determining the eigenvalues of an x 
matrix which can be explicitely calculated from the initial data for any discrete-time value 
n. However, one may ask the legitimate question what one has gained from this, as one 
has reduced the problem of solving one set of algebraic equations (given by the discrete 
equations of motion ( p.l4| )) to another problem of solving an algebraic equation (namely 
the factorization of a characteristic determinant). In the first case, nonetheless, we have 
to solve the system step by step at each iteration of the map, whence the issue of the 
complexity of the solution (i.e. the growth of the number of iterates) becomes imminent. 
In view of the integrability of the correspondence, one can at most have only polynomial 
growth, as Veselov has proven in [25|, and this is immediately clear from the exact solution. 



because the only growth comes from the permutations of eigenvalues. In fact, there are 
A^! permutations at each iteration of the map, but whatever be the permutation chosen at 
a certain step, it will end into one of the A^! possible permutations at the following step: 
the branches will "cross" each other with no divergence. 

i) Rational Case 

This is the limit that both periods tend to infinity, i.e. 2uji — > oo , 2012 — > ioo, in which 
case we can make the substitutions 

a{x) X , C,{x) \ , p{x) \ . 



X x^ 



Then, the Lax matrices take the form 



L^ = i/i*/i + Lo , M^ = -h% + Mo, (4.1) 

K K 
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where h denotes the (column-) vector with entries hi, the (row-) vector */i being its trans- 
posed, and in which 



TV 



TV 



»j=l A + g» - .^.^^ A + q^-qj 

From the form of the Lax matrices (4.2), we can then derive the following relations 



(4.2) 



AMo + QMo - MoQ = h% , 
XLq + QLq - LqQ = h% , 



(4.3a) 
(4.3b) 



where we have set 



TV 

Q — ^ ^ Qi^ii ■ 
i=l 



(4.4) 



On the other hand, from the Lax equation (2.7) and inserting ([4.1[), we obtain the relations 



LoMo = MqLq , 
Lfjh — Moh = —ph , 



(4.5a) 
(4.5b) 
(4.5c) 



together with the conservation law % ■ h = ■ h. Consequently, we can put 

Mo = UoUq^ , Lo = UoAUo^ , 



(4.6) 



where Uq is an invertible N x N matrix, and where the matrix A is constant, A = A , as 
a consequence of ( |4.5a ). We assume that A can be chosen to be a diagonal matrix, but 
this is not essential for what follows. Then, introducing 



'h-Uo , 



(4.7) 



we obtain from (^^) and 



(pi + A) ■ r = r , *s • (p/ + A) = *s 



(4.8) 



in which / denotes the N x N unit matrix, as well as from (4.3a) and ( 4.3b| ) 

X + V -V = ¥^s , AA + [F,A1 = r*s , 



(4.9) 
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together with the conservation law r^s = r^s . Eliminating the dyadic r *s from eq. ( [4.S| ) 
by making use of ( [4.8D , we find the linear equation 

V = {pl + AyWipI + A) - -f^ , (4.10) 



which can be immediately solved to give 

pX 



V{n) = (p/ + A)-" 



V{0) - n- 



{A + pir . (4.11) 



pI + A_ 

Thus, the solution of the discrete RS model in the rational case is given by the following 
statement : the position coordinates {qi{n)} of the particles at the discrete time n, evolving 
under the discrete equations of motion l{2. 1^ ) (in the rational limit cr(x) ^-^ x), are given 
by the eigenvalues of the matrix: 

Q{0)-np\{pI + Lom-^ , (4.12) 

where p can be chosen equal to unity without loss of generality. 

The choice of p only affects a scaling of the variables hi in (2.1), which enter the initial 
value of the matrix Lo(0)- 

ii) Hyperbolic (Trigonometric) Case 

In the hyperbolic limit 2uji oo , 2uj2 = ^Tri , in which case we can make the substitutions 

a{x) sinh(x) , ({x) — > coth(x) , p{x) — > sinh~^(x) , 

(after an appropriate gauge transformation of the Lax matrix with a diagonal matrix), 
leading to 

$K(a;) coth(K) + coth(2;) . (4.13) 

In this case we obtain for the Lax matrices ( p.l[ ) 

= Lo + (cothK + 7) , = Mo + {coth K + j)h% , (4.14) 

where 7 can be chosen at our convenience. In this case the reduced Lax matrices Lq and 
Mq are given by 

JV 

-^0 = 51 (coth(gi - qj + A) - 7) Cij , (4.15a) 

N 

Mo= J2 ^i^j (coth(^i - qj + A) - 7) Cij . (4.15b) 
12 



For them, as a consequence of the 'sphtting off' of the terms with the spectral parameter 



cothK, we get again the system of equations ( [4.5D . From eqs. (4-15) we derive that 

e2^e2Q(Lo + 7/1 - (Lo + -fh%)e'^Q = e^^e^^h^h + h^he^^ , (4.16a) 
g2Ag2Q(^^ + 7/1 - (Mo + 7/1 %e^^ = e^^e^^h % + h , (4.16b) 

in which Q, h and */i are given as before. We again make the identifications (^^), and 
introduce V,r,s as in ([4.7|), for which we subsequently derive the relations 

g2Ag2y (^A + (7 - l)r *s) = (a + (7 + l)r *s) e^^ , (4.17a) 
g2Ag2y + _ i)f tg) = + + i)f . (4.17b) 

Choosing now 7 = 1, and using (|4.8| ) we can again eliminate the dyadic r *s and explicitely 
solve for exp(2y). We obtain 

= {p + A)-^e^^{pe-^^ + A) , (4.18) 

which can be immediately integrated to yield 

g2y{n) = (A + p/)-"e2^(°)(A + pe'^^lT , (4.19) 

the V{Q) and ?7(0) as well as A being determined from the initial data (5(0) and Q{V). 
Thus, rewriting ( [4.19|) we obtain the following result: the exponentials of the position 
coordinates {e^'^»(")} of the particles at the discrete time n, evolving under the discrete 
equations of motion l\2. 14 ) (in the hyperbolic limit cr{x) 1— > sinh(x) J, are given by the 
eigenvalues of the matrix: 

(Lo(0) + p/)~"e2QW(Lo(0) + pe-^^ir , (4.20) 

where p can be chosen equal to unity without loss of generality. 



The choice of p only affects a scaling of the variables hi in (2^), which can be accounted 
for in the calculation of the hi(0) from eq. (2.13). We note again that at each discrete-time 
value the positions of the particles is uniquely determined up to a permutation. 

We finally remark that the trigonometric limit of the elliptic functions 2cji = ^vr, 
2lo2 ioo is integrated along similar lines after doing the replacements 

a{x) — > sinx , C{x) cot x . 

We shall omit the details. 
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5 Connection with Discrete Soliton Systems 



In 0, H, the discrete CM model was obtained from pole-solutions of a lattice version of 
the Kadomtsev-Petviashvili (KP) equation. The idea that integrable particle models of 
CM type are connected with soliton equations goes back already to the late seventies, cf. 



[P6| and also |27, but had never been applied to discrete soliton equations. In [10|, 
but maybe a bit more transparently in a connection between the relativistic particle 
model and soliton solutions of nonlinear integrable PDE's was established. It is a natural 
question to ask whether such a connection also exists on the discrete level. 

In order to establish such a connection, let us investigate more closely the trigonometric 
solution of the previous section. After a gauge transformation with a diagonal matrix of 



the form P 



(A + p/)(A + e^'^^pl) 



1/2 



we can transform e^^ = P "■WP"' into 



in which pi = (P" • = (*s • e^^P~")i , and uj = e^'^. The dependence on the discrete 
(time-) variable is given by 



-1 _\ 1/2 



leading to the equations for W 



^ Xi+u '^p , 

Pi = — Pi > (5-2) 



AW -ujWA = -2p^p , (5.3a) 
PW - ujWP = -2p^p , (5.3b) 

We can then introduce the characteristic polynomial 

N 

riO = Y[i^ + e^'n = dem + W) , (5.4) 

i=i 

and show that t(^) plays the role of the tau-function of a discrete soliton system. In fact, 
W is the kernel of the integral operator in the soliton sector that stands at the basis of 
the so-called direct linearization approach, cf. e.g. p0|]-p2[. 



To derive the relevant equations directly from the resolvent of the matrix W, we 
proceed as follows. First, using eq. (|5.3b| ), we can perform the following simple calculation 

m 



, ^, det(iujCI + W)-\CI + W 
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det 



uj-^ det (l + {uCI + W)-\A + pi)-' "* 



leading to 



viioC) = l+'p- {lo^I + W)-\K + pi)-' ■ p = u;""^^ . (5.5) 



which will turn out to be one of the variables governed by the soliton system. In order to 
derive discrete soliton equations for r(^) or v, we introduce the A^-component vectors 

ujiO ^ + W)-'A^ ■ p , 'ujiO ^ V • A'i^I + W)-' , (5.6a) 

as well as the scalar variables 

Ui,i^) = 2*p • A^^I + W)-'A' ■ p , (5.6b) 



G Z). Making use of the relations ( ^.3D , we can derive the following set of recursive 
relations between the different vectors Uj(^) and *Uj(^) 

ojP ■ UjiO = puj{uj-'C) + Uj+i{uj-'0 + , (5.7a) 

%i^) ■ P = Lo'u,+i{tuC) +p'Mu;0 - 'uoiu^UojiO ■ (5.7b) 



Using the definition ( |5.6b ) as well as ( |5.2D , we can derive also 



pUijiO + Ui+uiO = pUijiiuO + ujUi j+iiioO - Ui,o{uj^)UojiO ■ (5.7c) 

At this point we need to reflect a moment on the role of the discrete-time shift. The shift 
^ijiO ^ijiO is 'labeled' by the variable p, which can be identified as the (reciprocal) 
of the lattice parameter. Let us assume that the variables Uij{S,) do not depend on only 
one discrete time-variable n, but on a number of them, say n^, (a = 1, 2, 3, . . .), and that 
the corresponding discrete flows are compatible, each associated with its own parameter 
p, i.e. pi,P2,P3, .... This means that ?7 is a function of these independent time- variables, 
i.e. Uij{(,) = Uij{S,;ni,n2, . . .). In that case we have for instance that Uij{S,;ni,n2, ■ ■ ■) = 
ni + 1, n2 + 1, . . .) = Uij{^; ni, n2, . . .), etc. , the ^ corresponding to the translation 
in the second discrete variable. Of course, this puts extra compatibility conditions on 
the eigenvalues e^'^^ of the kernel W given in ( ^.Ij ) of the soliton solutions. In fact, the 
particle-coordinates qj can be shown to obey some two-dimensional lattice equations in 
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terms of the shifts ~ and ^, and the compatibiUty of such equations were investigated in 
for the non-relativistic situation. Here, we restrict ourselves to obtaining from ( |5.7(j ) 
closed-form nonlinear partial difference equations for special elements or combinations of 
elements of Uij, i.e. partial difference versions of the well-known soliton systems. These 
difference systems were investigated in detail in a number of earlier papers, (cf. |^3| for 
a review, and references therein), on the basis of the system of equations ( |5.7aD -( ^.7(j ). 
Thus here we will only present the results, referring to those papers for their derivation. 

In fact, for the special element Uni,n2,nz 

711,712,^3), we can derive 

the partial difference equation 

(Pl P2 ~\~ ^ni+l,n2,n3 'Uni,n2+l,i3 'Uni+l,n2+l,n3 

) + cycl. perm. = , (5.8) 



which is a lattice version of the Kadomtsev-Petviashvili (KP) equation, cf. [30, 31|. In 
[^, § pole solutions of this equation were investigated associated with a discretization of 
the non-relativistic Calogero-Moser model. The relation to the r-function Tn-^^n2,n2. = 
r(a;~*^"^"'""'2+"-3)^. j^^^ ^2, 713) is given by 

I / \ ''"rai+l,n2+l,n3''"rai, 12,713 /r n\ 

PI-P2 + Ifni,n2+l,n3 " n„^+l,„2,n3 = (Pl - P2) ■ ■ , (5.9) 

'Tni+l,n2 yU^'Tni ,n2+l,7i3 

from which one can derive the bilinear equation 

(Pl -P2)-rni+l,n2+l,n3('^~^0Tni,n2,n3+l(0 + cycl. perm. = 0, (5.10) 



which was first presented in (sj], cf. also |35]. Eq. ( 5.10| ) is related also to the following 



version of the lattice (modified) KP equation, p5|], 

/ \ I f^71l ,712 + 1,713 + 1 ^71i + l,n2 + l,7l3 \ / \ / ^71.l,712 + l,n3 + l 'y71i+l,7l2,n3 + l \ 

[Pl -P2) = [Pi -Ps) , 

y ^711 ,n2, 713+1 'i^ni+l,7i2,n3 / y "^711,712+1,713 ^711+1,712,713 / 

(5.11) 

in terms of the variable ^m, 712,713 = f (w"^"^"''"'^'''"'^-'.^; 7ii, 712, 713) , with v being defined in 
(|5^ , (taking pi=p). 



Finally, we mention that special reductions of the above equations arising by imposing 
additional symmetries on the soliton solutions, will lead to the lattice soliton systems of 
Gel'fand-Dikii type that were introduced in [^2|. In fact, taking the parameter uj equal to 
a root of unity, we can derive additional constraints on the system of equations given by 
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(|5^) -(|5?7c|). Thus, we find in the particular case of = —1 for u the following lattice 



version of the Korteweg-de Vries equation 

{Pl -P2 + Wni,n2+1 - Uni+l,n2){Pl + P2 + Um,n2 " ""m+l.na+l) = Pi - pI , (5-12) 



cf. [32 1, which is related to the following bilinear equation in terms of the r-function 

{Pi + + {P1 P2)Tni — l,n2'^ni+l,n2+l — '^PlTn\,n2'^n\,n2+l • 

(5.13) 

It is this equation that corresponds in the continuum limit exactly to the special case 



considered in [29|. Similar equations can be derived for other values of to when |a;| = 1, 



starting from the results presented in |3 



6 Connection with Bethe Ansatz Equations 

We already remarked above that the equations of motion of the discrete Ruijsenaars- 



Schneider model, (2.14), resemble closely the form of Bethe Ansatz equations (BAE's) for 
certain integrable quantum models. In this section we will make this connection more 
precise. Let us first focuss on the hyperbolic limit of eqs. ( ^.14 ), and connect it to the 



BAE's for the XXZ spin-^ Heisenberg magnet, cf. e.g. |36]-[^ and references therein. 



After that, we will show that the general form of ( |2.14| ) in the elliptic case is connected to 
the generalized Bethe Ansatz (BA) for the Heisenberg XYZ model, proposed first in 
and developed further in |^^. We will not treat the correspondence in the rational (or 
XXX) case since it can be obtained by a simple limit from the hyperbolic case. 

Let us, thus, consider the quadratic i?-matrix algebra 

(u — V )T(i)(n)r(2)(i;) = r(2)(^;)T(i)(n)i?(i2) 

{u- v) , (6.1) 

which is one of the central relations in the quantum inverse scattering method (QISM), (cf. 



for an early, and |38, 41] for more recent reviews), and let us consider the i?-matrix 
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for the spin-i XXZ Heisenberg magnet, 



/ a \ 

6 c 
c 6 

\ a y 

in which 

a = sinh(ii + ?]) , 6 = sinhu , c = sinhr7. 
The quantum L-operator on each site k of the spin chain has the form 

sinh(ti + ^rjal) sinh(r/)(T^ 



(6.2) 



sinh(77)(T^ sinh(n — ^rja^) 



(6.3) 



where = ^(cr^ ± icj|) , cr^'^'^ being the Pauh matrices on site k, (i.e. o"^ = 1 (g) . . . (g) 
fj" . . . (g 1 with 0"" on the kth entry of the tensor product). Let us now construct the 
following monodromy matrix 

/ 



T{u) = Lm{u - 5m) ■..L2{u- 52)Li{u - Si 



(6.4) 



V 



A{u) B{u) 
C{u) D{u) 

where M is the length of the spin chain, and in which the 6k,k = 1, . . . , M are impurity 
parameters. Both the Lk{u) as well as T{u) obey the relation (|6.1| ). Let us now recall the 
standard algebraic Bethe Ansatz construction, |3^]-|^, i.e. there is a vacuum state VLm-, 
which is an eigenvector of the diagonal entries A{u), D{u) of the monodromy matrix, and 
which is annihilated by the operator C{u) 



3.5) 



As the monodromy matrix is a comultiplication of the L-operators along the sites of the 
chain, we have that 

M M 

- - (6.6) 



M M 

o-iu) = Y[ ak{u) , d{u) = Jl dk{u) , 

k=l k=l 

where the functions ak{u) and dk{u) are the eigenvalues of the vacuum for the diagonal 



entries of the corresponding L-operator on site k. The algebraic Bethe Ansatz amounts 
to the creation of an eigenstate of the form 



N 



,qN) 



Yl B{qj)nM 



(6.7) 
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of the trace of the monodromy matrix A(u) + D{u), 

{A{u) + D{u))^ = tiu)^ , (6.8) 
with t{u) being the corresponding eigenvalue. We have then the following proposition: 



Proposition: ^I' is an eigenfunction of the trace of the monodromy matrix iff the num- 
bers Qj satisfy the Bethe Ansatz equations 

^ sinh(,,-,. + ,) ^aM^ £=l,...,iV. (6.9) 
smh(gf£ -qk-V) d{qe) 



In the case of the monodromy matrix (6.4), the functions a{u), d{u) take the form 



M ^ M ^ 

aiu) = n sinh('u - 5k + -r]) , d{u) = J| sinh(n -6k- -rj) , (6.10) 

k=l k=l 

which in the special case that M = 2N lead to the BAE's 

AT . , , AT s'mli{qi -Qk) smh.{qi - Qk+X) 

^T7 u = 11 ^T7 r~T7 = u ' ^=l>--->^> (6-11) 

smh(g^ -qk- X) ^J-^ smh(g^ - gk) smh{qe - qk - X) 



i.e. eq. ( 2.14| ) in the hyperbolic limit ( (y{x) ^ sinh(x) ) with the identifications A 



r/, qk = 5k — ^r], qk = dk+N + ^Vi = 1,...,A^). Thus, the discrete RS system in 
the hyperbolic limit can be reinterpreted as the BAE's for the spin-^ XXZ Heisenberg 
chain with the number of spins (impurities) equal to twice the excitation number of the 
eigenstate. We should stress that within such an identification the impurities play the 
role of the particle coordinates at times n — 1 and n + 1 {q resp. q), which means 
that the integrable correspondence mixes the parameters of the quantum model (i.e. the 
impurities) and the parameters of a solution given by the Bethe wavefunction. 

To make the connection with the elliptic case we have to consider the QISM for the 
XYZ Heisenberg model that was treated in |39| for the spin-^ situation. For arbitrary 
spin one needs to consider the algebras introduced by Sklyanin in [^]. However, we are 
interested here in the inhomogeneous spin-i chain, i.e. including arbitrary impurities. 
Both extensions were treated in a recent paper by Takebe [^^. It was shown that the 
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u 



(6.12) 



generalized Bethe Ansatz of Takhtajan and Faddeev is applicable to this model. The 
corresponding i?-matrix is given by 

/ a d \ 
6 c 
c 6 
\ d a J 

in which 

a = sn{u + 2r]; k) , b = sn{u;k) , c = sn(27/; /c) , d = ksn{2r]; k)sn{u; k)sn{u + 2r]] k) . 
The quantum L-operator on each site k of the spin chain now has the form 



Lfc(u) 



wqI + w^al wial - iw2<7l 



(6.13) 



where 



W0 + W3 = sn(n + ri;k) , wq - = sn(ii - ri;k) , 
wi + W2 = sn(2T/; k) , — t(;2 = A;sn(2r/; A;)sn(ii + r/; A;)sn(u — 77; k) 



The monodromy matrix takes again the form ( |6.4| ). The generalization of the algebraic 
Bethe Ansatz to the model described in ( |6.12| ), ( |6.13 ) is known only for the special eigen- 
functions for which the number of excitations is half the number of spins in the chain. 
Of course, this constitutes only a special class within the total set of eigenfunctions, but 
no further results exist to date. Curiously, this is precisely the case in which we have a 
connection with the discrete RS model! 

To be more precise, the generalized BA, developed in []39| , |40[| , consists of constructing 
eigenfunctions of the trace of the monodromy matrix A{u) + D{u) after performing a 
gauge transformation of the form 

/ 



Tin) 



'{u) = Sm{u) T{u)Sm'{u) =: 



(6.14) 



^m.m' (^) ^m,m'(^) 
y Cm,m' {u) DjYi.m' (u) 

where Sm{u) is an appropriate scalar matrix. (Note that we deviate from the notations in 



9|, 4C] to be consistent with the notations used in earlier sections of the present paper.) 
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Instead of a single vacuum there is now a set of generating (vacuum) vectors i^M,n for 
which one has 

AM,n{u)i^M,n = a{u)^lM,n-l , DM,n{u)^M,n = d{u)Q,M,n+l , CM,n{u)^M,n = . 

(6.15) 

The generahzed Bethe eigenfunctions are hnear combinations of the vectors 



^n(9l, • • • ,qN) = -Bn+l,n-l(9l) • • • Bn+N,n-NiqN)^M,n-N , 



(6.16) 



of the form 



(6.17) 



where is a free parameter. As it was pointed out in |39], it is possible that the series 
( |6.17| ) is summable to zero for all 9 except for a finite number of values 9j . (The results of 



Baxter |36] show that among the 6j there is also the value 9 = 0.) Then for such values of 
9 one can produce the corresponding eigenfunctions, provided that the excitation numbers 
qi,...,qN obey the following BAE's (cf. [||]) 



V TT ^(^^ -qk + ^rj) a{qi - 6k + 2r/) 

^ 11 _ _ 9,.^ - 11 



tJi (^{qe -qk- 2??) 



k=l 



cr{qe - h) 



1, 



X = exp 



N 



2N 



UJ2 



, k=l 



k=l 



(6.18) 



(6.19) 



in the special case^ that M = 2N, and re-expressing the equations derived in |^0| in terms 



of cr-function. Eqs. ( 6.18 ) correspond to the equations ( 2.14| ) for the discrete RS model 
in the elliptic case by the following identifications: A = 27], qk = 6^, qk = ^k+N — 2??, 
(A; = 1, . . . , N), together with 

N 



P 

- = exp 
P 



Am9 - A— [B + B-2B 



k=l 



(6.20) 



The sum of the particle positions B plays the role of the boost generator of the underlying 
Poincare algebra, cf. [|l^. So, the generalized Bethe Ansatz equations correspond to the 
discrete-time system with the dynamics of p given in terms of the dynamics of qkS. The 



^ The relation M = 2N follows because of spin-i situation, while becoming A'^ = £M for the arbitrary 
spin i (cf. H). 
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results of the Section 3 can be generalized to incorporate this case. In fact, we need to 
add now the term 

A7ri9B - ^X(B - Bf , (6.21) 

to the Lagrangian ( p. lb ) in order to get the correct factor from the Euler-Lagrange equa- 
tions. If we now calculate the canonical momenta from ( 3.1b| ) we get an extra term 
-^A(5 - 1 

.2 : 



^^A(i? — B) on the right-hand side in ( p. 41 ), but this factor cancels again when we 



express the hj in terms of the canonical momenta, because of the p entering into the eqs. 

Thus, the final result Q remains unaltered^. Consequently, we still have the 
same integrals in terms of the canonical variables, and the Ruijsenaars variables do not 
change, hence they are in involution and we have the integrable correspondence. 

It is amazing that exactly in the restricted case where the generalized BA is known 
to apply, i.e. the case that the number of impurities is equal to twice the number of 
excitations, we have a connection with an integrable correspondence of RS type. This 
might also suggest that there exists a deeper relation between integrable quantum models 
solvable by the QISM on the one hand and classical integrable multiparticle models on the 
other hand. The connection provided here might yield an explanation for the existence 
conditions of the generalized BA of Faddeev and Takhtajan and of Baxter's solution of 
the eight-vertex model. 



7 Discussion 



The connection exhibited in the previous section between the BAE's for the XXZ and XYZ 
model with impurities and the discrete relativistic multi-particle model is very remarkable. 
A similar connection seems to exist also in the non-relativistic case studied in [|l|, |[ , which 
yields a link between the equations of motion of the discrete CM model and BAE's for 
the Gaudin model. In both cases, we thus obtain a new interpretation of the BAE's for a 
quantum solvable model, namely as a classical integrable dynamical model with discrete 
time. Such interpretations, together with the exact solutions of the discrete equations 
of motion that were derived in section 4, might lead to possible ways of 'solving' the 
■^This effect can be viewed as a consequence of the relativistic invariance of the modeL 
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BAE's, (which in principle constitutes a set of transcendental equations), without having 



to perform some thermodynamic limit, cf. |38]. This, however, requires a reformulation 
of the initial value problem for the discrete model rather as a boundary value problem on 
the discrete-time chain. Although it is not clear at this stage that such boundary value 
problems are explicitely solvable as it was the case in Section 4, we nevertheless conjecture 
that this intriguing dynamical interpretation of the Bethe Ansatz equations will shed new 
light on the solvability of the corresponding quantum models. Furthermore, it would be 
interesting to study generalizations of our results, for instance in the directions of finding 
discrete particle models associated with higher-spin models and higher-rank associated 
Lie algebras (nested Bethe Ansatz), as well as investigating in this light the quadratic 
algebras associated with boundary conditions for quantum integrable systems (reflection 
equation algebras), cf. |43|. The latter might lead to time-discretizations of BCn type of 



RS models, which were investigated in [14, 15]. 
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